Multiple roles of brassinosteroid signaling in vascular development.

Brassinosteroids (BRs) are plant steroid hormones that control growth and stress responses. In the context of development, BRs play diverse roles in controlling cell differentiation and tissue patterning. The vascular system, which is essential for transporting water and nutrients throughout the plant body, initially establishes a tissue pattern during primary development and then dramatically increases the number of vascular cells during secondary development. This complex developmental process is properly regulated by a network consisting of various hormonal signalling pathways. Genetic studies have revealed that mutants defective in BR biosynthesis or the BR signalling cascade exhibit a multifaceted vascular development phenotype. Furthermore, BR crosstalk with other plant hormones, including peptide hormones, coordinately regulates vascular development. Recently, the involvement of BR in vascular development, especially in xylem differentiation, has also been suggested in plant species other than the model plant Arabidopsis thaliana. In this review, we briefly summarize the recent findings on the roles of BR in primary and secondary vascular development in Arabidopsis and other species.

The evolution and functional divergence of FT-related genes in controlling flowering time in Brassica rapa ssp. rapa.

KEY MESSAGE: The BrrFT paralogues exhibit distinct expression patterns and play different roles in regulating flowering time, and BrrFT4 competes with BrrFT1 and BrrFT2 to interact with BrrFD proteins. Flowering time is an important agricultural trait for Brassica crops, and early bolting strongly affects the yield and quality of Brassica rapa ssp. rapa. Flowering Locus T paralogues play an important role in regulating flowering time. In this study, we identified FT-related genes in turnip by phylogenetic classification, and four BrrFT homoeologs that shared with high identities with BraFT genes were isolated. The different gene structures, promoter binding sites, and expression patterns observed indicated that these genes may play different roles in flowering time regulation. Further genetic and biochemical experiments showed that as for FT-like paralogues, BrrFT2 acted as the key floral inducer, and BrrFT1 seems to act as a mild 'florigen' protein. However, BrrFT4 acts as a floral repressor and antagonistically regulates flowering time by competing with BrrFT1 and BrrFT2 to bind BrrFD proteins. BrrFT3 may have experienced loss of function via base shift mutation. Our results revealed the potential roles of FT-related genes in flowering time regulation in turnip.

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f.

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

Limited associations between MHC diversity and reproductive success in a bird species with biparental care.

The selective pressure from pathogens on individuals can have direct consequences on reproduction. Genes from the major histocompatibility complex (MHC) are central to the vertebrate adaptive immune system and pathogen resistance. In species with biparental care, each sex has distinct reproductive roles and levels of investment, and due to a trade-off with immunity, one can expect different selective regimes acting upon the MHC of each parent. Here, we addressed whether couples combine each other's variation at MHC loci to increase their breeding success. Specifically, we used a 23-year dataset from a barn owl population (Tyto alba) to understand how MHC class Iα and IIß functional divergence and supertypes of each parent were associated with clutch size and fledging success. We did not detect associations between MHC diversity and supertypes with the clutch size or with the fledging success. In addition, to understand the relative contribution from the MHC of the genetic parents and the social parents, we analyzed the fledging success using only a cross-fostered dataset. We found several associations of weak-to-moderate effect sizes between the father's MHC and fledging success: (i) lower MHC-Iα divergence in the genetic father increases fledging success, which might improve paternal care during incubation, and (ii) one and two MHC-IIß DAB2 supertypes in the social father decrease and increase, respectively, fledging success, which may affect the paternal care after hatching. Furthermore, fledging success increased when both parents did not carry MHC-IIß DAB1 supertype 2, which could suggest conditional effects of this supertype. Although our study relied on a substantial dataset, we showed that the associations between MHC diversity and reproductive success remain scarce and of complex interpretation in the barn owl. Moreover, our results highlighted the need to incorporate more than one proxy of reproductive success and several MHC classes to capture more complex associations.

The effects of co-invasion by three Asteraceae invasive alien species on plant taxonomic and functional diversity in herbaceous ruderal communities in southern Jiangsu, China.

Invasive alien species can affect plant taxonomic and functional diversity. Multiple invasive alien species can co-invade the same plant community. However, the effects of such co-invasion on plant taxonomic and functional diversity are currently unclear. Our study aimed to estimate the effects of co-invasion by three Asteraceae invasive alien species (i.e., Conyza canadensis (L.) Cronquist, Conyza sumatrensis (S.F. Blake) Pruski and G. Sancho, and Solidago canadensis L.) on plant taxonomic and functional diversity in herbaceous ruderal communities in southern Jiangsu, China. The effects of these three invasive alien species under seven invasion combinations (including invasion by one invasive alien species, co-invasion by two invasive alien species, and co-invasion by these three invasive alien species) on plant taxonomic and functional diversity were investigated in a comparative field study of herbaceous ruderal communities. Niche differentiation mediated the functional divergence between these three invasive alien species and natives under all invasion combinations. These three invasive alien species significantly increased plant taxonomic diversity (especially plant diversity and richness) and plant functional diversity (especially Rao's quadratic entropies) under all invasion combinations. The relative abundance of invasive alien species was significantly positively associated with plant functional diversity (especially community-weighted mean trait values and Rao's quadratic entropy). The number of invasive alien species was significantly positively associated with plant taxonomic diversity (especially plant diversity and richness) and plant functional diversity (especially Rao's quadratic entropies). Thus, co-invasion by these three invasive alien species may synergistically increase plant taxonomic diversity (especially plant diversity and richness) and functional diversity (especially Rao's quadratic entropies).

The evolution and functional divergence of 10 Apolipoprotein D-like genes in Nilaparvata lugens.

Apolipoprotein D (ApoD), a member of the lipocalin superfamily of proteins, is involved in lipid transport and stress resistance. Whereas only a single copy of the ApoD gene is found in humans and some other vertebrates, there are typically several ApoD-like genes in insects. To date, there have been relatively few studies that have examined the evolution and functional differentiation of ApoD-like genes in insects, particularly hemi-metabolous insects. In this study, we identified 10 ApoD-like genes (NlApoD1-10) with distinct spatiotemporal expression patterns in Nilaparvata lugens (BPH), which is an important pest of rice. NlApoD1-10 were found to be distributed on 3 chromosomes in a tandem array of NlApoD1/2, NlApoD3-5, and NlApoD7/8, and show sequence and gene structural divergence in the coding regions, indicating that multiple gene duplication events occurred during evolution. Phylogenetic analysis revealed that NlApoD1-10 can be clustered into 5 clades, with NlApoD3-5 and NlApoD7/8 potentially evolving exclusively in the Delphacidae family. Functional screening using an RNA interference approach revealed that only NlApoD2 was essential for BPH development and survival, whereas NlApoD4/5 are highly expressed in testes, and might play roles in reproduction. Moreover, stress response analysis revealed that NlApoD3-5/9, NlApoD3-5, and NlApoD9 were up-regulated after treatment with lipopolysaccharide, H2 O2 , and ultraviolet-C, respectively, indicating their potential roles in stress resistance.

Functional Divergence and Origin of the Vertebrate Praja Family.

The Praja family is an E3 ubiquitin ligase, promoting polyubiquitination and subsequent degradation of substrates. It comprises two paralogs, praja1 and praja2. Prior research suggests these paralogs have undergone functional divergence, with examples, such as their distinct roles in neurite outgrowth. However, the specific evolutionary trajectories of each paralog remain largely unexplored preventing mechanistic understanding of functional differences between paralogs. Here, we investigated the phylogeny and divergence of the vertebrate Praja family through molecular evolutionary analysis. Phylogenetic examination of the vertebrate praja revealed that praja1 and praja2 originated from the common ancestor of placentals via gene duplication, with praja1 evolving at twice the rate of praja2 shortly after the duplication. Moreover, a unique evolutionary trajectory for praja1 relative to other vertebrate Praja was indicated, as evidenced by principal component analysis on GC content, codon usage frequency, and amino acid composition. Subsequent motif/domain comparison revealed conserved N terminus and C terminus in praja1 and praja2, together with praja1-specific motifs, including nuclear localization signal and Ala-Gly-Ser repeats. The nuclear localization signal was demonstrated to be functional in human neuroblastoma SH-SY5Y cells using deletion mutant, while praja2 was exclusively expressed in the nucleus. These discoveries contribute to a more comprehensive understanding of the Praja family's phylogeny and suggest a functional divergence between praja1 and praja2. Specifically, the shift of praja1 into the nucleus implies the degradation of novel substrates located in the nucleus as an evolutionary consequence.

Structure and functional divergence of PIP peptide family revealed by functional studies on PIP1 and PIP2 in Arabidopsis thaliana.

PAMP-induced secreted peptide (PIP), one of the small post-translationally modified peptides (PTMPs), plays a crucial role in plant development and stress tolerance. However, little is known about functional divergence among this peptide family. Here, we studied the evolution of the PIP family in 23 plant species (10 monocotyledons and 13 dicotyledons from 7 families) and their functional divergence in Arabidopsis. A total of 128 putative PIP precursors were identified and classified into two subfamilies through phylogenetic analysis. Functional studies on AtPIP1 which represents Clade I family and AtPIP2 which represents Clade II family have shown that AtPIP2 displayed stronger immunity induction activity but weaker root growth inhibition than AtPIP1 in Arabidopsis. Transcriptome analysis of Arabidopsis seedlings treated with AtPIP1 and AtPIP2 showed that differential genes for both polypeptides were significantly enriched in similar plant defense pathways. However, Co-expression and Protein-protein interaction (PPI) analysis showed that the functions of AtprePIP2 co-expressed genes were more enriched in plant defense pathways than AtprePIP1. Molecular docking results show that AtPIP1 binds to RLK7 receptor with a more stable free energy and less binding area than AtPIP2, while hydrogen bond transfer occurs at the SGP motif position. The above results suggest that the PIP family have undergone functional divergence during evolution. Collectively, this work illustrates the relationship between PIP structure and function using Arabidopsis PIP as an example, and provides new insights into the current understanding between growth inhibition and immune responses which may be correlated but not fully coupled.

Molecular Evolution and Protein Structure Variation of Dkk Family.

Dkks have inhibitory effects on the Wnt signaling pathway, which is involved in the development of skin and its appendages and the regulation of hair growth. The nucleotide sequences were compared and analyzed to further investigate the relationship between the structure and function of the Dkk gene family and vertebrate epidermal hair. The analysis of the molecular evolution of the Dkk family revealed that the evolution rate of the genes changed significantly after speciation, with the Aves and Reptilia branches showing accelerated evolution. Additionally, positive selection was observed at specific sites. The tertiary structure of the protein was also predicted. The analysis of the functional divergence of the Dkk family revealed that the functional divergence coefficient of each gene was greater than 0, with most of the functional divergence sites were located in the Cys-2 domain and a few in the Cys-1 domain. This suggests that the amino acid and functional divergence sites may play a role in regulating the binding of the Dkk family to LRP5/6, and thus affect the inhibition of Wnt signaling, leading to different functions of Dkk1, Dkk2, and Dkk4 in the development of skin hair follicles. In addition, the Dkk families of Aves and Reptilia may have undergone adaptive evolution and functional divergence.

Parallel evolution of picobirnaviruses from distinct ancestral origins.

IMPORTANCE: Picobirnaviruses (PBVs) are highly heterogeneous viruses encoding a capsid and RdRp. Detected in a wide variety of animals with and without disease, their association with gastrointestinal and respiratory infections, and consequently their public health importance, has rightly been questioned. Determining the "true" host of Picobirnavirus lies at the center of this debate, as evidence exists for them having both vertebrate and prokaryotic origins. Using integrated and time-stamped phylogenetic approaches, we show they are contemporaneous viruses descending from two different ancestors: avian Reovirus and fungal Partitivirus. The fungal PBV-R2 species emerged with a single segment (RdRp) until it acquired a capsid from vertebrate PBV-R1 and PBV-R3 species. Protein and RNA folding analyses revealed how the former came to resemble the latter over time. Thus, parallel evolution from disparate hosts has driven the adaptation and genetic diversification of the Picobirnaviridae family.

Evolution and functional role prediction of the CYP6DE and CYP6DJ subfamilies in Dendroctonus (Curculionidae: Scolytinae) bark beetles.

Dendroctonus-bark beetles are natural components and key ecological agents of coniferous forests. They spend most of their lives under the bark, where they are exposed to highly toxic terpenes present in the oleoresin. Cytochrome P450 (CYP) is a multigene family involved in the detoxification of these compounds. It has been demonstrated that CYP6DE and CYP6DJ subfamilies hydroxylate monoterpenes, whose derivatives can act as pheromone synergist compounds or be pheromones themselves in these insects. Given the diversity and functional role of CYPs, we investigated whether these cytochromes have retained their function throughout the evolution of these insects. To test this hypothesis, we performed a Bayesian phylogenetic analysis to determine phylogenetic subgroups of cytochromes in these subfamilies. Subgroups were mapped and reconciled with the Dendroctonus phylogeny. Molecular docking analyses were performed with the cytochromes of each subgroup and enantiomers of α-pinene and ß-pinene, (+)-3-carene, ß-myrcene and R-(+)-limonene. In addition, functional divergence analysis was performed to identify critical amino acid sites that influence changes in catalytic site conformation and/or protein folding. Three and two phylogenetic subgroups were recovered for the CYP6DE and CYP6DJ subfamilies, respectively. Mapping and reconciliation analysis showed different gain and loss patterns for cytochromes of each subgroup. Functional predictions indicated that the cytochromes analyzed are able to hydroxylate all monoterpenes; however, they showed preferential affinities to different monoterpenes. Functional divergence analyses indicated that the CYP6DE subfamily has experimented type I and II divergence, whereas the CYP6DJ subfamily has evolved under strong functional constraints. Results suggest cytochromes of the CYP6DE subfamily evolve to reinforce their detoxifying capacity hydroxylating mainly α- and ß-pinene to (+) and (-)-trans-verbenol, being the negative enantiomer used as a pheromone by several Dendroctonus species; whereas cytochromes of the CYP6DJ subfamily appear to retain their original function related to the detoxification of these compounds.

Enhanced plasticity and reproductive fitness of floral and seed traits facilitate non-native species spread in mountain ecosystems.

Floral and seed traits, their relationships, and responses to abiotic constraints are considered the key determinants of the invasion success of non-native plant species. However, studies evaluating the pattern of floral and seed traits of non-native species in mountain ecosystems are lacking. In this study, we determined (a) whether the floral and seed traits of native and non-native species show similarity or dissimilarity across elevations in mountains, and (b) whether the non-native species follow different allometric patterns compared with native species. Functional variations between native and non-native species were assessed through floral and seed traits: flower count, flower display area, flower mass, specific flower area, seed count, and seed mass across an elevational gradient. Permanent plots (20 × 20 m) were laid at each 100 m elevation rise from 2000 to 4000 m a.s.l. for sampling of herbaceous plant species. The mean values of floral and seed traits such as flower display area, specific flower area, and seed count were significantly higher for non-native species compared to native species. A significant difference in trait values (flower display area, flower mass, seed count, and seed mass) between non-native species and native species was observed along the elevational gradient, except for flower count and specific flower area. The bivariate relationship revealed non-native species to exhibit a stronger relationship between flower display area ∼ flower mass, and flower display area ∼ seed mass traits than the native species. Non-native species showed enhanced reproductive ability under varying environmental conditions along an elevational gradient in mountain ecosystems. Greater flower display area and seed mass at lower elevations and a stronger overall trait-trait relationship among non-native species implied resource investment in pollinator visualization, flower mass, and seed quality over seed quantity. The study concludes that enhanced plasticity and reproductive fitness of floral and seed traits would consequently aid non-native species to adapt, become invasive, and displace native species in mountain ecosystems if the climatic barriers acting on non-native species are reduced with climate change.

Expression profiles in knock-down transgenic plants of high and low diversified duplicate genes in Arabidopsis thaliana.

Duplicated genes show various degrees of functional diversification in plants. We previously identified 1,052 pairs of high diversified duplicates (HDDs) and 600 pairs of low diversified duplicates (LDDs) in Arabidopsis thaliana. Single knock-down of HDDs induced abnormal phenotypic changes because the other gene copy could not compensate for the knock-down effect, while single knock-down of LDDs did not induce abnormal phenotypic changes because of functional compensation by the copy gene. Here, focusing on one pair each of HDDs and LDDs, we performed transcriptome analyses in single-knock-down transgenic plants. The numbers of differentially expressed genes in single-knock-down transgenic plants were not different between HDDs and LDDs. Thus, functional compensation inferred by transcriptomics was similar between HDDs and LDDs. However, the trend of differentially expressed genes was similar in the pair of LDDs, while expression profiles were dissimilar in the pair of HDDs. This result indicates that a pair of LDDs tends to share similar functions but a pair of HDDs tends to have undergone functional divergence. Taking these findings together, as the reason for no phenotypic changes in single knock-down of LDDs but phenotypic changes in double knock-down of LDDs, we concluded that phenotypic changes of LDDs were induced by decreasing gene dosage.

Components of functional diversity revisited: A new classification and its theoretical and practical implications.

Functional diversity is regarded as a key concept for understanding the link between ecosystem function and biodiversity. The different and ecologically well-defined aspects of the concept are reflected by the so-called functional components, for example, functional richness and divergence. Many authors proposed that components be distinguished according to the multivariate technique on which they rely, but more recent studies suggest that several multivariate techniques, providing different functional representations (such as dendrograms and ordinations) of the community can in fact express the same functional component. Here, we review the relevant literature and find that (1) general ecological acceptance of the field is hampered by ambiguous terminology and (2) our understanding of the role of multivariate techniques in defining components is unclear. To address these issues, we provide new definitions for the three basic functional diversity components namely functional richness, functional divergence and functional regularity. In addition, we present a classification of presence-/absence-based approaches suitable for quantifying these components. We focus exclusively on the binary case for its relative simplicity. We find illogical, as well as logical but unused combinations of components and representations; and reveal that components can be quantified almost independently from the functional representation of the community. Finally, theoretical and practical implications of the new classification are discussed.

Host functional traits as the nexus for multilevel infection patterns.

Understanding pathogen transmission and infection patterns at multiple biological scales is a central issue in disease ecology and evolution. Here, we suggest that functional traits of host species readily affect infection patterns of species, communities, and landscapes, and thus serve as a linkage for multilevel studies of infectious disease.

Functional traits explain waterbirds' host status, subtype richness, and community-level infection risk for avian influenza.

Species functional traits can influence pathogen transmission processes, and consequently affect species' host status, pathogen diversity, and community-level infection risk. We here investigated, for 143 European waterbird species, effects of functional traits on host status and pathogen diversity (subtype richness) for avian influenza virus at species level. We then explored the association between functional diversity and HPAI H5Nx occurrence at the community level for 2016/17 and 2021/22 epidemics in Europe. We found that both host status and subtype richness were shaped by several traits, such as diet guild and dispersal ability, and that the community-weighted means of these traits were also correlated with community-level risk of H5Nx occurrence. Moreover, functional divergence was negatively associated with H5Nx occurrence, indicating that functional diversity can reduce infection risk. Our findings highlight the value of integrating trait-based ecology into the framework of diversity-disease relationship, and provide new insights for HPAI prediction and prevention.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat.

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Review: WRKY transcription factors: Understanding the functional divergence.

WRKY transcription factors (TFs) play crucial roles in the growth and development of plants and their response to environmental changes. WRKY TFs have been detected in sequenced plant genomes. The functions and regulatory networks of many WRKY TFs, especially from Arabidopsis thaliana (AtWRKY TFs), have been revealed, and the origin of WRKY TFs in plants is clear. Nonetheless, the relationship between WRKY TFs function and classification is unclear. Furthermore, the functional divergence of homologous WRKY TFs in plants is unclear. In this review, WRKY TFs were explored based on WRKY-related literature published from 1994 to 2022. WRKY TFs were identified in 234 species at the genome and transcriptome levels. The biological functions of ∼ 71 % of AtWRKY TFs were uncovered. Although functional divergence occurred in homologous WRKY TFs, different WRKY TF groups had no preferential function.

Climate change threats to the global functional diversity of freshwater fish.

Climate change impacts on freshwater ecosystems and freshwater biodiversity show strong spatial variability, highlighting the importance of a global perspective. While previous studies on biodiversity mostly focused on species richness, functional diversity, which is a better predictor of ecosystem functioning, has received much less attention. This study aims to comprehensively assess climate change threats to the functional diversity of freshwater fish across the world, considering three complementary metrics-functional richness, evenness and divergence. We built on existing spatially explicit projections of geographical ranges for 11,425 riverine fish species as affected by changes in streamflow and water temperature extremes at four warming levels (1.5°C, 2.0°C, 3.2°C and 4.5°C). To estimate functional diversity, we considered the following four continuous, morphological and physiological traits: relative head length, relative body depth, trophic level and relative growth rate. Together, these traits cover five ecological functions. We treated missing trait values in two different ways: we either removed species with missing trait values or imputed them. Depending on the warming level, 6%-25% of the locations globally face a complete loss of functional diversity when assuming no dispersal (6%-17% when assuming maximal dispersal), with hotspots in the Amazon and Paraná River basins. The three facets of functional diversity do not always follow the same pattern. Sometimes, functional richness is not yet affected despite species loss, while functional evenness and divergence are already reducing. Other times, functional richness reduces, while functional evenness and/or divergence increase instead. The contrasting patterns of the three facets of functional diversity show their complementarity among each other and their added value compared to species richness. With increasing climate change, impacts on freshwater communities accelerate, making early mitigation critically important.

Gene family expansion and functional diversification of chitinase and chitin synthase genes in Atlantic salmon (Salmo salar).

Chitin is one of the most abundant polysaccharides in nature, forming important structures in insects, crustaceans, and fungal cell walls. Vertebrates on the other hand are generally considered "nonchitinous" organisms, despite having highly conserved chitin metabolism-associated genes. Recent work has revealed that the largest group of vertebrates, the teleosts, have the potential to both synthesize and degrade endogenous chitin. Yet, little is known about the genes and proteins responsible for these dynamic processes. Here, we used comparative genomics, transcriptomics, and chromatin accessibility data to characterize the repertoire, evolution, and regulation of genes involved in chitin metabolism in teleosts, with a particular focus on Atlantic salmon. Reconstruction of gene family phylogenies provides evidence for an expansion of teleost and salmonid chitinase and chitin synthase genes after multiple whole-genome duplications. Analyses of multi-tissue gene expression data demonstrated a strong bias of gastrointestinal tract expression for chitin metabolism genes, but with different spatial and temporal tissue specificities. Finally, we integrated transcriptomes from a developmental time series of the gastrointestinal tract with chromatin accessibility data to identify putative transcription factors responsible for regulating chitin metabolism gene expression (CDX1 and CDX2) as well as tissue-specific divergence in the regulation of gene duplicates (FOXJ2). The findings presented here support the hypothesis that chitin metabolism genes in teleosts play a role in developing and maintaining a chitin-based barrier in the teleost gut and provide a basis for further investigations into the molecular basis of this barrier.